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Abstract 

We study in this work the properties of the Q m f network which is constructed from an 
anisotropic partition of the square, the multifractal tiling. This tiling is build using a single 
parameter p, in the limit of p —> 1 the tiling degenerates into the square lattice that is asso- 
ciated with a regular network. The Q m f network is a space-filling network with the following 
characteristics: it shows a power-law distribution of connectivity for k > 7 and it has an high 
clustering coefficient when compared with a random network associated. In addition the Q m f 
network satisfy the relation N oc £ d f where £ is a typical length of the network (the average 
minimal distance) and N the network size. We call df the fractal dimension of the network. In 
the limit case of p — > 1 we have df — > 2. 

PACS numbers: 89.75.Da, 89.75.Hc, 61.43.Hv 
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I - INTRODUCTION 



The last years have seen an increasing interest in network studies in physics 0, 0. 
Despite graph theory have been a research topic of mathematics and science of compu- 
tation, the physicists have driven their attention to networks that show distribution of 
connectivity, P(k) : following a power law and that present small world effect: £ oc InN 
for £ a typical network distance and TV the size of the network. A recent article |3| points 
the difficulty to put together these two aspects in a common broad scale-free framework, 
it means, a fractal paradigm. In fact, a real fractal should have TV oc £ df for df the fractal 
dimension, and not a logarithm dependence between TV and £. In this work we present 
a network that have the following characteristics: it is fractal, TV oc £ d f , it is scale-free 
P(k) oc k~ 7 for a large range of fc, and it is complex in the sense of having large clustering 
coefficient. The scaling relation N oc £ d is trivially fulfilled for regular lattices where d is 
an integer, the topologic dimension of the space. The complex fractal network we explore 
in this paper has a non integer df. 

Recently networks embedded in metric spaces have been investigated in the literature 
because of the large applications of networks that effectively occupy a volume in 3- 
dimensions. In addition, motived by microprocessors design space-filling networks [5] have 
been studied. We cite these new trends in networks because the distinguished network 
we analyze in this paper is embedded in a metric space and furthermore it is space filling. 
In fact our network has a geometrical inspiration, it comes from a partition of the square. 
Indeed this network is originated from a singular tiling that has the additional property 



of being a multifractal partition of the^square 6|. 

0, B uM 



We follow the previous literature 0, 0, 0, Iq^IioI] and call this object the multifractal 
tiling, Q m f. The Q m f tiling was developed in the context of modeling transport and 
percolation in heterogeneous porous media. In fact, a broad set of irregular and het- 
erogeneous systems are model in the literature using a multifractal approach. We cite 
systems in geology HQ , atmospheric science an d chemistry [lj]. Oil reservoirs 
are complex anisotropic structures whose treatment have been challenged science and 
technique because of its non trivial geometry. Inspired in the description of oil reservoirs 
it was developed the Q m f tiling. In reference 0] an exhaustive study of the percola- 
tion threshold of Q m f was performed, in [8] a random version of Q m f was created, in 
some of its percolation critical exponents were found, and in a numerical study of its 
coordination number is done. 

In this work we explore the network properties of the Q m fi it means its topology, 
the study of the connections (neighborhood) among the cells of the tiling. The paper is 
organized as follows. In section 2 we describe in some detail the process of construction 
of the Q m f tiling and show its more important properties. In section 3 we show the 
main results concerning the properties of the network: the distribution of connectivity, 
the scaling of the the average minimal distance between two sites and the clustering 
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FIG. 1: The first five steps in the construction in the multifractal tiling Q m f for (s, r) = (1, 3). 

coefficient. Finally in section 4 we present our conclusions, discuss the main implications 
of our results and compare the properties of the Q m f network with other networks in the 
literature. 

II - THE MULTIFRACTAL OBJECT 

The multifractal tiling is a peculiar partition of the square. It is interesting to think 
about the Q m f tiling in contrast to the square lattice. The square lattice can be con- 
structed using the following algorithm: take a square and cut it symmetrically with 
horizontal and vertical lines. This procedure produces four square cells. Repeat this 
procedure n times inside each new block and you have finally a square lattice with 2 n 
cells. The Q m f object is generated in a similar way as the square lattice above described, 
but instead of using a symmetric partition we perform horizontal and vertical sections 
following a given ratio. 

In Fig. [I] we exemplify the five initial steps of the construction of the multifractal for 
the parameter p = 1/3, or (s, r) = (1,3). In [I] (a) we show, n — 0, the initial square 
that by convenience we assume of size L — 1. In (b), n — 1, a vertical cut is performed 
and two rectangles are formed. We call this a (5, r) = (1,3) object because the square 
is divided in 4 parts such that 1 part stays at one side and 3 parts at the other side. In 
(c), n — 2, two horizontal lines are drawn using the same section rate as before. At this 
level the initial square generates four rectangular blocks. Using as the area unit a square 
of size e = l/(s + r), the largest block has area r 2 , there are two blocks of areas rs and 
the smallest block has area s 2 . In (d) and (e), n — 3 and n — 4, resp ectively, the same 
procedure is repeated inside the initial four blocks. In reference |l5| it is explored the 
group of eight possibilities of cutting, a square lattice with a given ratio. In this work, as 
in other papers about Q m f 0, Q, Q Jill it is followed the recipe of Fig. [IJ 

We remark that the number of blocks at step n is 2 n . These blocks do not have all 
the same area, we call the subsets of blocks of same area by a fc-set. It is easy to check 
that the block area distribution follows a binomial distribution and the number of /c-sets 



3 




FIG. 2: A Q m f multifractal tiling and a Q m f network, for p = 1/2 and n = 6. In (a) the 
original tiling and in (b) the corresponding network. 



is k — n + 1. This fact implies that the Qm/ has the remarkable property: in the limit of 
n —> oc the area of its forming blocks follows a multifractal distribution p . The spectrum 
of fractal distributions comes from a box counting reasoning: 

log N(X) 

D x lim - - (1) 

e-o log (l/ e ) v ; 

for N(X) the number of unitary cells of size length e that cover the set of blocks of a 
given area X. Once the initial square is partitioned n times, the size of the unitary cell 
is e = l/(s + r) n . For each fc-set the total area of blocks (using e area units) is done by: 

N k = C k n s k r {n ~ k \ (2) 

where C k is the binomial coefficient that express the number of elements fc-type, and 
s k r( n_fc ) is the area of each element of this set. We put together all these elements to 
have the fractal dimension of each fc-set: 

logN k \og{C k s k r^) 
D k = lim - — — = hm — . (3) 

e^0log(l/e) n-^oo log (s + r) n KJ 

This distribution show a concave shape with a maximum at k = pn. The case r = s = 1 is 
degenerated. In this situation the subsets of the lattice are composed uniquely by square 
cells of the same area. Therefore the tiling is formed by a single subset of dimension 2. 
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FIG. 3: The cumulative degree distribution for several values of p and N — 2 . For k > 7 
the distributions approaches a power-law form. The exponent 7 controlling the decay of the 
distributions decreases in the limit of p — » 0, the full set of exponents is show in Table HI 



In Fig. El (a) it is shown an example of Q m f construction for (s, r) = (1, 2) and n = 6. 
In Fig. [21 (b) we build a network corresponding to this tiling. The nodes of the network 
are the 2 n blocks of the Q m f and the vertices are established according to a neighborhood 
criterium. These last figures offer a glimpse of the metric heterogeneity and the topology 
of the multifractal. In the next section we explore in detail the network properties of this 
class of objects. 



Ill - RESULTS 

We start the analysis of the properties of the Q m f network discussing its distribution 
of connectivity, P{k). In Fig. El we show the cumulative sum of P{k) versus k for several 
values of p as indicated in the figure. The option for the cumulative sum instead of 
P(k) itself is due to the strong fluctuation of the data. Fig. El confirms the results of a 
previous work For low fc, typically k < 7, the curve suggests an exponential behavior 
and above this threshold the network depicts a scale- free behavior. The values of the 
exponents 7 of the power-law P(k) oc k~ J are indicated in Table [J for several p, trivially 
the values of 7 are estimated from the slopes of the curve of the cumulative sum decreased 
by one. The exponent 7 goes to an asymptotic limit for large N |H|. We observe that in 
the limit of p — ► 1 the Q m f tiling gets more symmetric and at p — 1 the Q m f degenerates 
into the square lattice. For the regular square lattice P(k) is a delta of Dirac centered at 
k = 4 and the cumulative sum a step function. Fig. El corroborates this idea, the skewed 
curves in the p — > 1 limit anticipate the phase transition at p — 1. 

We explore the distance characterization of the Q m f network in Fig. El where we 
display the behavior of the average shortest distance for all couple of distinct vertices 
of the network, versus network size, N. The simulation is performed for some values 
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FIG. 4: In (a) it is displayed the average distance £ as a function of N for several values of p. 



of p as indicated in the figure. For all p examined the data show a power-law behavior 
N oc £ d f. We find 2 < df < 4, the full set of df is shown in Table [U Two limit cases are 
interesting. The limit p — > 1, which corresponds to the square lattice, has df — > 2 as it 
is expected in a bidimensional space. The opposit limit p —> 0, which is associated with 
very anisotropic structures, shows large df. We cannot affirm that 4 is an asymptotic 
threshold, further numerical investigation should test this hypothesis. We remark that 
the Qmf network does not follow a small world relationship £ oc In TV that is common to 
most of power-law and random like networks. 

An analysis of the clustering coefficient, C, versus network size, AT, is shown in Fig. El 
The general view of this figure points to a stable behavior of C in the limit of large N. 
The dispersion of C among p is not large, the numerics show C = 0.37 ± 0.01. Smaller 
values of p, however, show a significant larger C. The discussion about C is intriguing 
once we compare the numerical values of C with the clustering coefficient of a random 
network associated to the Q m f network. An associated random network is defined as a 
network with the same N and < k > of the original network (we do not compare our 
results with a random network with a same P(k) because such random network would 
alterate the space filling characteristics that we are interested in) . For a random network 
C =< k > /N : in the case of our network: < k > is a constant number smaller than 6 
and N a number that can grow without limit. As a consequence the associated random 
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FIG. 5: In the clustering coefficient C as function of N for the same set of values of p of the 
previous figure. 



network has C — > in the limit TV — > oc. Therefore the Q m / network has a C that is 
infinitely larger than the clustering coefficient of the associated random network. Because 
of the high C and the power-law behavior of distribution of connectivity we call the Q m f 
network a complex network. 

In Table U it is shown some parameters related to the Q m f network for several values 
of p. Most of these data was already commented in the text. We focus now on the 
average connectivity < k > of the network. For all p studied we have < k >~ 5.43 which 
characterizes a sparse network. This is not surprising, since there is a result in topology 
that shows that for two dimensions the average coordination number of a tiling cannot 
exceed 6. The average connectivity confirms Fig. El where we can see that the majority of 
vertices are situated in the range: 4 < k < 6. Otherwise we note that most of interesting 
results concerning the distribution of connectivity, in special the power-law behavior, are 
satisfied only in the range k > 7. Indeed, the Q m f network, as most of complex networks, 
also have hubs that determine the distinguished characteristics of the network. 
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TABLE I: The average quantities: clustering coefficient, C, minimal distance, £, and connectiv- 
ity, < k >; the slope 7 of the distribution of connectivity and the fractal dimension df. The 
data corresponding to the average parameters are estimated for N = 2 16 . 



p 


(3,4) 


(2,3) 


(1,2) 


(1,3) 


(1,4) 


(1,8) 


c 


0.3603 


0.3610 


0.3703 


0.3735 


0.3769 


0.3777 


£ 


122.02 


115.26 


96.43 


75.12 


57.12 


41.70 


< k> 


5.4357 


5.4357 


5.4338 


5.4275 


5.4357 


5.4357 


7 


16.6 


11.1 


6.4 


4.2 


3.4 


2.9 


df 


2.09 


2.14 


2.25 


2.52 


2.94 


3.79 



IV - CONCLUSION 

In this work we explore some properties of a space filling network that come from a 
multifractal partition of the square lattice, the Q m f network. An analysis of the distribu- 
tion of the connectivity, P{k) ) assures that the Q m f network shows a power-law tail that 
is more accentuated as increases the anisotropy of the underlying Q m f tiling. Roughly 
the power-law tail of P(k) starts at k ~ 7. We remark that there is no regular lattice in 
2 dimensions with k > 6 and typical Voronoi latices have an exponential small number 
of vertices in this range. In addition the Q m f network has a clustering coefficient that 
approaches a constant, C = 0.37 ± 0.01, that does not depend on N. This fact is in 
contrast to random networks that (for a constant < k >) have C oc AT -1 . Because the 
value of C is much larger than the value of C of the associated random network we call 
the Q m f network a complex network. 

The most interesting aspect of the Q m f network concerns its fractal behavior. For 
the average minimal path £ we observe that N oc £ d f for the fractal dimension, df. The 
simulations show that 2 < df < 4. The lower limit correspond to the case p — 1 where the 
multifractal tiling degenerates into the square lattice. In this situation the slope 7 (from 
the power law P(k) oc k~ J ) increases dramatically, this situation corresponds to the P(k) 
of the square lattice that is of the form of a Delta of Dirac. The opposite limit, p — 0, 
corresponding to very anisotropic tilings, presents comparatively small values of 7. 

We point that, diversly from [3], the Q m f network shows an actual fractal behavior A^ oc 
£ d f that is obtained without any renormalization artefact. In the reference [3] a ingenious 
procedure is used to calculate two fractal dimensions ds and df. The first dimension 
depends on a suitable embedding in a metric space and a box counting methodology. 
The second is based on network distance and a mass (number of vertices) inside a given 
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radios. In our case, we have calculated df in the standard way, the box-counting, however, 
depends on the methodology we use to make the embedding of the Q m f object. The 
simplest embedding is the Q m f lattice itself that is a 2-dimensional object and as a result 
ds = 2. Note that the multifractal property of Q m f appears when we consider the subsets 
of blocks of same area, if we disregard the area set a bidimensional tilling assumes the 
trivial topologic dimension, df = 2. 

The Q m f tiling is indeed a remarkable mathematical object, from a metric perspective 
it is a multifractal: it is formed by a denumerable quantity of sets of different areas each 
one with a given fractal dimension. In a topologic perspective the connections among the 
vertices (the cells of the tiling) form a fractal network. We remark that regular networks 
(generated from lattices for instance) and the Bethe tree satisfy the criterium TV oc £ d , but 
these structures are regular. For the best knowledgement of the authors the Q m f network 
is the only case of a true fractal, scale- free and with high clustering coefficient. 
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